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Abstract 

We argue that Landau- Fermi liquids do not have any gravity duals in the purely 
classical limit. We employ the logarithmic behavior of entanglement entropy to 
characterize the existence of Fermi surfaces. By imposing the null energy condition, 
we show that the specific heat always behaves anomalously. We also present a 
classical gravity dual which has the expected behavior of the entanglement entropy 
and specific heat for non-Fermi liquids. 
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1 Introduction and Summary 



The AdS/CFT correspondence pQ provides us a very useful tool to study various proper- 
ties of strongly coupled systems. It is clear that among many possible states in condensed 
matter physics, systems with Fermi surfaces are remarkably important. Recently, there 
have been exciting progresses on the holographic realizations of Fermi surfaces using the 
AdS/CFT such as the ones using the Reissner-Nordstrom black holes pHH] and the star- 
like solutions, called electron stars [HHH]. Refer also to [9H12] for further developments and 
other constructions. The former construction leads to the non-Fermi liquids, while the 
latter to the standard Landau-Fermi liquids. The non-Fermi liquids have anomalous be- 
haviors, for example, with respect to the specific heat and conductivity when we compare 
them with the Landau-Fermi liquids. 

An outstanding advantage of AdS/CFT is that we can analyze strongly coupled sys- 
tems in a classical gravity by ignoring all quantum corrections. However, there have not 
been found any holographic constructions of Fermi surfaces where basic quantities such 
as the free energy, the entropy and the specific heat are all calculated in the classical 
gravity and agree with the expected standard behaviors of Fermi surfaces. For example, 
a charged black hole whose near horizon includes AdS2 has the unusual property of the 
non-vanishing entropy even at zero temperature. On the other hand, in the electron star 
solutions, which approach to the Lifshitz geometry [13] in the IR limit, the specific heat 
and entropy are not large enough to be calculated in the classical gravity. The purpose 
of this paper is to study if it is possible to construct a holographic dual of Fermi surfaces 
in a purely classical gravity. If such a gravity dual exists, we expect that the number 
of Fermi surfaces is of order N 2 , taking the Fermi energy to be order one. Here we are 
considering a system where fermions belong to the adjoint representation of the SU(N) 
gauge group. The fractionalized fermi liquid [9] is a promising candidate of such 0(N 2 ) 
fermi surfaces. 

For this purpose, we need to make the definition of Fermi surfaces clear. In this paper, 
we employ the entanglement entropy [HHlTj to characterize the presence of Fermi surface. 
We define systems with Fermi surfaces by requiring that their entanglement entropies show 
the logarithmic violation of the area law. This has been shown explicitly for free fermion 
theories [IS] and recent progresses in condensed matter physics support that this property 
is true even for non- Fermi liquids [19,20j. 

By using the holographic entanglement entropy [2TH23] . we can rewrite this require- 
ment as the IR behavior of the metric in the gravity dual. We mainly focus on 2+1 
dimensional systems dual to asymptotic AdS4 geometries. Then we require the null en- 
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ergy condition, which is expected to be satisfied by any physically sensible backgrounds 
with well-defined holographic duals. This leads to additional constraints on the metric. 
Eventually we find that to satisfy these two requirements the specific heat C should be- 
have like C oc T a with the constraint a < 2/3 at low temperature T for 2+1 dimensional 
systems. Thus the standard Fermi liquid (a = 1) is not allowed. In this way, we find that 
under our classical gravity dual assumption, we can construct only non- Fermi liquids. We 
also give an example of an effective gravity action which has such a gravity solution. We 
will leave the embeddings of our gravity backgrounds in string theory as a future problem. 

This paper is organized as follows: In section two, we discuss the behavior of entangle- 
ment entropy in the presence of Fermi surfaces. In section three, we study the constraints 
on the gravity dual metric for the Fermi surfaces from the viewpoint of the entanglement 
entropy. In section four, we further impose the null energy condition and see how it con- 
strains the behavior of the specific heat. In section five, we present an effective gravity 
model which has the expected behavior of the entanglement entropy and specific heat for 
non-Fermi liquids. 

2 Fermi Surface and Entanglement Entropy 

Consider the entanglement entropy in a d+1 dimensional quantum field theory (QFT) on 
R 1 ' . The coordinate of R 1 ' is denoted by (t, X\,X2, Xd)- The entanglement entropy Sa 
associated with a subsystem A is defined by the von- Neumann entropy Sa = — Trp^ log pA- 
Pa is the reduced density matrix defined in terms of the total density matrix p by pa = 
Tr bPi where the subsystem B is the complement of A. 

We choose the subsystem A as a strip with the width /. This is explicitly expressed as 

A = {(x!,x 2 , ■ ■ -,x d )\ - - < xi < < x 2 ,x 3 , • • -,x d < L}. (2.1) 
In a d + 1 dimensional conformal field theory, the entanglement entropy Sa behaves 

like 

^ = 7^-^, (2.2) 

where e is the UV cut off (e — > in the UV limit); 7 and a are numerical constants which 
represent the degrees of freedom [T6l[T7] . This has been confirmed in the holographic 
calculation [2Tir2~2~] . The leading divergent term is proportional to the area of the boundary 
of A and this is called the area law P3J[T5]. This is true for all QFTs with UV fixed points. 
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2.1 Logarithmic Violation of Area Law 

On the other hand, for ad+1 dimensional QFT with a Fermi surface, we expect that 
the finite part of Sa is substantially modified when the size I of the subsystem A is large 
enough 

S * = ^^TT + riL d ~ l k d f 1 \og(lk F ) + 0(1°), (2.3) 

where kp is the Fermi momentum and rj is a positive numerical constant. 

Before we give a field theoretic explanation of the behavior (12.31) below, we would like 
to point out that (12. 3ft is essentially the same as the logarithmic violation of area law in 
femionic theories, which is well known in the condensed matter literature. 

In the papers [18] (see also the review [H]), the violation of area law in free fermions 
on lattices has been shown in any dimension d. They argue that the leading term in the 
entropy looks like 

S A ~L d ~ 1 \ogL+---, (2.4) 

where we assume that the subsystem A is given by a finite size (but large) region (e.g. a 
d dimensional ball) which spans order L in any directions. Their calculation is based on 
a discretized lattice setup, and there they assume that the Fermi surface has a finite area. 
In our continuum limit of the quantum field theories, this corresponds to the case where 
the Fermi energy is the same order of the UV cut off scale i.e. kp ~ e _1 . Indeed, in this 
limit, the second term in H 2 . 3 1) gets dominant and leads to the violation of the area law 
(12 .4p by setting I ~ L. However, in this paper what we are interested in is the case where 
kp is finite while e is taken to be zero. Therefore we have the expression (12. 3 p instead of 
flUD- 

Even in the presence of strong interactions this kind of the violation of area law has 
been expected. The recent calculation of Sa for a critical spin liquid state j2()] can be 
regarded as an evidence of such a log term even in a strongly coupled fermion theory with 
a Fermi surface. Also such a violation of area law has also been argued in |19j . 

2.2 Free Fermion Calculation 

To understand ( 12.31) . we would like to consider a calculation of S A for a free fermion in 
(d+l) dimension We compactify on a torus with a large radius L. The 

momentum in these directions are quantized as 

ki=-, (i = l,2,.-.,d). (2.5) 
4 Similar estimations have already been done in Q1J1H3]. See also [21] for an explicit calculation. 
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We assume there exists a Fermi surface with the Fermi momentum Jzf and the ground 
state is describes as 

l*> = II & Ln 2 ,-,nJO>. (2.6) 
|fc|<fc F 

Since the theory is free, we can treat sectors with different quantum numbers n = 
(ri2,ri3,- • -jTid) as decoupled independent sectors. Each sector with a fixed number of 
n can be regarded as a 1 + 1 dimensional free massive fermion theory with the Fermi 



surface at k\ = y k\ — k 2 . The mass is given by m = \k\. The ground state in the 1 + 1 
dimensional theory for a fixed n looks like 

i**>= n d°>> ( 2 - ? ) 

ki<y/k%-\k\* 

and the total ground state is given by 

l*> = IIl*»>- ( 2 - 8 ) 

n 

Notice that when kp < \k\ there is no Fermi surface in the reduced two dimensional 
fermion theory. 

Since the density matrix is factorized into sectors with various n, we can express Sa 
as the summation over each sector: 

4 d+1) = E4 1+1) (-)- (2-9) 

n 

We are interested in the region 

/ » ^, (2.10) 

kp 



so that the entanglement entropy can probe the IR physics. When \k\ > kp, there is no 
Fermi surface and S A 1+1 \n) can be found from the known result in two dimensional CFT. 
If we define S A to be the entanglement entropy for a massive two dimensional QFT 
without any Fermi surface, then it behaves like 

S% FT2 ~ logZ/e (Im « 1), 

S<f T2 ~ -logme (lm»l). (2.11) 



Note also that we omitted the coefficient, which is proportional to the central charge |16j . 
Since I » ^=-, when \k\ > kp we find 

S% +1) (n) ~-log(\k\e). (2.12) 
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In this way we we obtain the following estimation 

c^+D ^ Ld -t r {dk) d - x s F A (\k\A + L d - 1 [ {dk) d ^s F A (\%\, Je F - \kA , 

Jk F <\k\<e-^ ^ J J\k\<k F \ J 

~ -L d - 1 [ (dk) d - 1 log(e\k\) + L d - 1 [ {dkY^S^ (\k\, Jk 2 F -\k\A , (2.13) 

ifc F <|fc|<e- 1 J\k\<k F V / 

where S A (m, fi) is the entanglement entropy for a free massive fermion in 1 + 1 dimension 
with the mass m and the Fermi momentum /i. The first term in ( 12 . 1 3[) comes from the 
modes \k\ > kp can be simplify estimated as 

-L d - X [ (dfc^logCelfcl)- (-Y V--, (2.14) 

which just represents the standard law. 

To proceed further we need to know S A for non-zero values of /i. The dispersion 
relation in the massive two dimensional fermion obtained by the dimensional reduction 
looks like 



Ek = \ \k\ 2 + [\^k F -\k\ 2 + 5ki 



h 2 — \h\ 2 
rv j? I n/ I 



^ k F + ~ ; Sh. (2.15) 

Kp 

Therefore near the Fermi surface (i.e. Ski << kp), the fermion behaves like a massless 
fermion. Therefore we can estimate the remained part of S A as follows 



L d - 1 [ (dk) d - l S F A (\k\, Jk 2 F -\k\ 2 ) 
J\k\<k F V / 



~ L d - 1 / {dk) d - l \ogl/t. (2.16) 

J\k\<k F 

This is because the low energy excitation around the Fermi surface becomes relativistic 
as in fl2~T5l) . 

Therefore, by combining (12. 14ft and (12 . 16p . we find that the total behavior of S A looks 
like _ 

S %+V~(^\ 1 + L ^k d F 1 \og(lk F ) + --: (2.17) 

Notice that the logarithmic divergent terms O(loge) are obviously canceled between (12. 14j) 
and (I2.16p . In this way, we managed to derive the formula ( 12. 3p . 
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So far we assumed the free fermions. In order to compare the results with those in 
AdS/CFT, we need to see how this argument changes due to the strong interactions. 
However, there are evidences which argue that the logarithmic dependence ( 12.31) does not 
change. One of them is the recent calculations done for the interacting Fermi surfaces in 
spin liquids [20J . Another is the fact that though in non- Fermi liquids we often encounter 
non-relativistic dynamical exponents for the dispersion relation around the Fermi surface, 
the logarithmic behavior (12.1 ip does not change (only its coefficient changes). Therefore, 
in this paper, we define systems with Fermi surfaces by requiring that their entanglement 
entropies show the logarithmic violation of the area law, as already mentioned. 



2.3 Entanglement Entropy for 2D Free Massless Dirac Fermions 

It is useful to look at an analytical expression of the entanglement entropy in two dimen- 
sional CFT at finite chemical potential. Below we will calculate Sa for a two dimensional 
free massless Dirac fermion on a circle at finite temperature T — 1//3 and at finite chemi- 
cal potential fi. We would like to ask the readers to refer to [25] for details where a similar 
analysis has been done when \i = 0. 

We consider a torus whose lengths in the space and the Euclidean time direction 
are given by (3 and 1, respectively. The size of subsystem A is denoted by L such that 
< L < 1. The thermal partition function is expressed by using the theta functions 



nh 



Tre -/3 J H'-27r/3 At 7V 



(2.18) 



where H and N are Hamiltonian and the fermion number. 

The two point function of twisted operators a k in the orbifolded free fermion theory 
on a torus is given by 



(<7 fc (L)<7_ fc (0)> 



2nr](i/3) 3 


4A fc 


0i{L\iP) 





9 3 (i/3fi\i/3) 



(2.19) 



where A fc = is the conformal dimension of the fc-th twisted operator a k . By applying 
the replica trick, we obtain 



Sa 







dN 



pogTt[(pA)"]] 







N=l 



dN 



JV-l 
2 



log n (a k (L)a_ k (0)} 



JV-l 
2 



(2.20) 



N=l 
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In the end, we find the high temperature expansion 



Sa(L,P,h) = -log 



— sinh ^ ~^ 



7T6 



m=l 



1 _ e 2 ^ e '^f )(1 - e -^f e -™f ) 



2tt# -2ttI 



(l-e- 2 ^)2 



/ 1 V / # COth ( # ) - 1 

+ 2^fc^(^ V^V I cos(2tt^), 



1=1 



I 



sinh 



7rA 



where e(-> 0) is the UV cut off. 

We find basic properties as follows 



(2.21) 



S A (L, (3, li) = S A (L, P,l-n) = S A (L, $,n + l), 



(2.22) 



which is explained by the particle-hole exchange symmetry and the periodicity of energy 
levels peculiar to the two dimensional free fermion theory. 
We can also confirm that the difference 



S A (L = l-5,f3,fi)-S A (L = S,(3,Li) 



7T 

3^ 



2£ 



-1^ 



nlL 



coth 



(*)- 



sinh ( f 



in the limit 5 — > actually coincides with the thermal entropy 

^(/3, / i) = -/3 2 ^(r 1 iog^). 



COs(27T/x/), 

(2.23) 
(2.24) 



In the low temperature expansion entanglement entropy is given by (for < /i < 1) 



S A (L,f3,fi) = -log 



— sin(7rL) 



Tie 



m=l 



(1 _|_ e -27T/3(m-l/2))2 



(-1)' /TrZLcot(TrLZ) - 1\ 
+ 2 ^ ^—^ „:„u V-/3f\ ) cosh(27r/i/3/) 



sinh (7T/3Z) 



(2.25) 



Thus at zero temperature, we simply find 



S a (L,oo,li) = - log 



sin(7rL) 



(2.26) 



we confirm that the logarithmic behavior still exists even if we turn on the chemical 
potential as we assumed in the previous calculations. 
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3 Possible Holographic Duals with Fermi Surfaces 



We are interested in holographic duals of strongly coupled 2 + 1 dimensional systems with 
Fermi surfaces. Thus we consider 3 + 1 dimensional gravity backgrounds. We consider 
the following metric 

ds 2 = ^ (-f(z)dt 2 + g(z)dz 2 + dx 2 + dy 2 ) . (3.27) 

This is the most general metric when we require the translational and rotational symmetry 
in x and y directions^ Notice that the metric component gt z can be eliminated by 
a coordinate transformation. We require that its dual theory has a UV fixed point and 
therefore (I3.27P should be asymptotically AdS 4 . Since z = corresponds to the boundary, 
we require 

/(0) = g(0) = 1. (3.28) 



3.1 Holographic Entanglement Entropy 

Now we holographically calculate the entanglement entropy in these gravity backgrounds. 
The holographic entanglement entropy [2TH23] is given by 

S A = ^>, (3.29) 

where 7^4 is the codimension two minimal area surface which coincides with OA at the 
boundary z = 0. 

We are interested in the strip shape subsystem A (12. ip i.e. 

A={(x,y)\-t<x<-, 0<y<L}. (3.30) 

We can specify the minimal area surface 7^ by the surface x = x(z) in (I3.27p . 
The area of this surface can be found as 



— ^g(z) + x'{z)\ (3.31) 



where z* is the turning point, where x' gets divergent; e is the UV cut off. The variational 
principle for x(z) leads to 

x'(z) 
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const. , (3.32) 



Here we did not assume any breaking of symmetry as we are interested in realizing a fermi liquid in the 
absence of any external fields. However, it is also possible that the logarithmic behavior of entanglement 
entropy occurs in the presence of external fields. A good example will be the magnetic quantum critical 
point found in an asymptotically AdSs geometry [26], as pointed out in [27] 
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which can be solved as 



The width I is related to by 



1=2 1 d4\ fz^- ^ 



In the end, the area (13.311) is expressed as follows 



dz / g{z) 



Area( 7A ) = 2R 2 L / — (3.35) 



^ 2 Vl-7 



To obtain explicit results, let us assume that the function g(z) scales as 

g(z) ~ (2 » 2 F ), 

~ 1 (z « (3.36) 

with a certain scale zp. The parameter Zj? is dual to the length scale where the geometry 
starts modified and shows peculiar IR behaviors. We assume n > 1 below. We are 
interested in the very IR limit z* » zp and would like to see the I dependence of the 
finite part of Sa- 

Then (13.341) is estimated as 

2 f z * z 2 ~^ n z n ^ 

I ~ — — / dz — ===== ~ c n — - — , (3.37) 

Z * Z F J a/1--2 4 /^ Zp 

where c n is a positive numerical constant. 

On the other hand, the area is evaluated as 

2R 2 L 2R 2 L f z * , z n ~ 2 



Area(7 A ) ^ 1 — / dz- . . — . 

e z F J ,Jl - z t/ z l 

2R 2 L , R 2 LzT 1 



+ d n (3.38) 

t Zp 

where d n is a positive numerical constant. 

In the end, the behavior of Sa is obtained as follows (k n is a positive constant) 

n-l 

R 2 L R 2 L ( iy* , 

S a = —jaT + fc «?r- • — +■■■, ( 3 - 39 ) 



2G$e G nZf \z f 
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where the omitted term represents the subleading term in the limit z* >> zp or equally 
/ >> zp. The leading divergent term agrees with the area law and this is expected because 
our background is asymptotically AdS. Note that the maximum increasing rate of Sa as a 
function of I is linear, which corresponds to n — > oo limit. This is consistent with the fact 
that the maximal allowed entropy is always given by the log of the dimensional of Hilbert 
space of A, which is clearly proportional to /. This limit is realized in the charge AdS 
black holes whose the near horizon geometry AdS2xR 2 . In this case, the spacetime has a 
horizon with non-vanishing entropy and the extensive behavior of Sa is indeed expected. 

On the other hand, for the purpose of a gravity dual of a Fermi surfaces, we need to 
set n = 1. In this case, we indeed find the following behavior (notice that (I3.37P remains 
the same, while we need to modify (13.381) ) 

Sa = + h^- log (-) + 0(f), (3.40) 

2G { ffe G N z F \z F J 

which agrees with (12. 3p . Therefore the existence of Fermi surfaces requires the behavior 




(z » z F ), 



~ 1 (z« z F ). (3.41) 

The parameter zp is now interpreted as the scale of the Fermi surface ~ k F . More 
precisely, since we are expecting the existence of many fermi surfaces, zp is regarded as 
the average of their fermi levels. This is our basic guiding principle to search a holographic 
dual of Fermi surfaces. It is also straightforward to extend the above argument to higher 
dimensions. 



3.2 Circular Subsystems 

So far we analytically confirmed the logarithmic dependence on the subsystem size I only 
for the strip subsystems. Here we would like to study the case where the subsystem A is 
given by a circular disk with radius I. In conformal field theories, we obtain 

S A = 7 • ~ - S, (3.42) 

where 7 and 5 are constants [T7J|2TJ[28]. In the presence of the Fermi surfaces, we again 
expect 

S A = l-- + f)-lkp\oglkF + ---. (3.43) 
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Figure 1: The shape of the function r(z) for = 1. 

In the holographic calculations with the metric (I3.27p . we introduce a polar coordinate 
such that dx 2 + dy 2 = dr 2 + r 2 d9 2 . The subsystem A of the CFT corresponds to r < I at 
z = in the gravity background. The surface 7^ is specified by r = r(z) in (13.271) at a 
fixed time t. The holographic entanglement entropy is given by minimizing 

Sa = ^ f* ^ r (s) v^iy+^W. (3.44) 

The boundary condition of r(z) is r(0) = / and r(z m ) = 0. At the turning point z = z*, 
r'(z) gets divergent. The equation of motion for r(z) is given by 

d ( r(z)r'(z) 
Z \z 2 ^g(z) + r'(z) 2 

The behavior near z = z* is completely fixed by the equation of motion ( 13. 45ft and given 
in the form 

r(z) ~ (z, - zf' 2 (r + ri(z„ - z) + • • •) . (3.46) 
We are working for the choice 

g(z) = V / l+^/4, (3.47) 

which is asymptotically AdS± and leads to the IR geometry ( 13.411) . Since we need to 
perform the numerical analysis below, we will simplify set = I. 

The explicit form of r(z) is shown in FigJTJ We plotted the finite part of the entan- 
glement entropy in Figj2]for the choice Zp — 1. Its behavior is indeed well approximated 
by 

S f A in ~ 0.50 -/log/ - 0.60-/. (3.48) 




y/g{z) + r'{zf 



(3.45) 
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Figure 2: The finite part of Sa (called S^ n ) as a function I. 

By looking at the behavior of the coefficients in (I3.48p . we can numerically confirm the 
dependence of Zf oc kp 1 as in (13.431) . 

4 Null Energy Condition and Specific Heat 

To have a physically sensible gravity theory, we need to require an appropriate energy 
condition. In the presence of negative cosmological constants, we usually impose the null 
energy condition: 

T^N»N U > 0, (4.49) 

where denotes any null vector. For example, we employ this condition to prove the 
holographic c-theorem [29J. 

We find from the Einstein equation that the energy stress tensor is simply given 

by 

Tpv = R„ u - -g^R. (4.50) 

We can choose as 

N l = — r =, N z = -y==, N x = sin9, (4.51) 

where 9 is an arbitrary constant. Then 

rp N , NV zg{z)f'{zf + f{z){zf{z)g'{z) + g(z)(4f\z) - 2zf"{z)) . a 

= AzJWgJzf Sm{d) 



zf(z)g(z) 
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and so the null energy condition (I4.49p is satisfied if and only if 

g(z)f'(z) + f(z)g'(z) < 0, 

zg(z)f(z) 2 + f(z)(zf'(z)g'(z) + g(z)(4f'(z) - 2zf"(z)) < 0. (4.53) 
Let us focus on the IR geometry. If we assume the behavior 

f{z) ex z~ 2m , g(z) cx z 2n , (4.54) 
then the conditions (14.531) lead to 

m > n. (4.55) 
4.1 Constraints on Specific Heat 

Since the Landau-Fermi liquid has the peculiar property that the specific heat is linearly 
proportional to the temperature, it is helpful to investigate the specific heat. At finite 
temperature, the IR geometry (I4.54p is modified into a black hole solution of the form: 



ds 2 = R 2 (-z- 2 ( m+1 >h(z)dt 2 + z 2 ^~h{z)- x dz 2 + dx2 + d y 2 ^j 



(4.56) 



where we set R = 1. We define the location of the horizon to be z = zu. For — << 1, the 
metric should approach ( 14.541) and therefore we have h(z) ~ 1 and h(z) ~ 1. On the other 
hand, in the near horizon region z ~ zh, we generically have the regular non-extremal 
horizon and thus they behave like 

h{z) ~ h(z) ~ fHJlf, (4.57) 

up to constant factors. 

Thus we can estimate the temperature as 

T oc z H m ~ n -\ (4.58) 

by requiring that the Euclidean geometry at z = zh is smooth. Since the thermal entropy 
S is given by the area of horizon we find 

S oc oc V 2 -T^+r, (4.59) 

Z H 

where V 2 is the volume of the dual CFT 3 . 
The specific heat is found as 

C=^ = T^ocT^n. (4.60) 
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The existence of Fermi surfaces requires n = 1 as explained in ( 13.41)) . In this case, the 
null energy condition f !4.55|) leads tc§ 



m + n + 1 3 

Since the standard Fermi liquid has the linear specific heat C oc T, classical gravity duals 
cannot have the Landau- Fermi liquids. Instead, our specific heat (14.601) suggests that they 
correspond to non-Fermi liquidgj. 

If the Fermi surface is govern by a two dimensional theory with the dispersion relation 
E ~ k z , then we expect 

CcxSocT^, (4.62) 
which agrees with the analysis in [32]. The null energy condition requires 

z > | (4.63) 

In summary, we found that the Landau-Fermi liquids cannot holographically be real- 
ized in any purely classical gravity. However, non-Fermi liquids are not excluded as long 
as the bound ( I4.63P is satisfied (in 2 + 1 dimension). For such a classical gravity dual, we 
expect the presence of order N 2 Fermi surfaces, taking the Fermi energy to be order one. 

It might be useful to note that the electron stars [&HE] are not included in our definition 
of Fermi surfaces with classical gravity duals. Actually, in an electron star, since there is 
a smaller number (<< N 2 ) of Fermi surfaces [7j, the entanglement entropy and specific 
heat contributed by them only appear as quantum corrections to the classical gravity 
calculations (see [8] for the analysis of the specific heat). On the other hand, the free 
energy is computable in classical gravity because the chemical potential is taken to be 
very large as 0(\/N) in the electron star. The same is also true for the IR gapped 
models |12|. 



5 Effective Gravity Model 

Finally, we would like to give an effective gravity model which has the logarithmic behavior 
( 12.31) of the holographic entanglement entropy, peculiar to systems with Fermi surfaces. 

6 In AdS5 gravity duals, we find C oc T m +™+ 1 and the null energy condition requires m > n. The 
Fermi surface (or equally the logarithmic entanglement entropy) corresponds to n = 2. Therefore we find 
the behavior C oc T a with the constraint a < | for the 3 + 1 dimensional fermi surfaces. 

7 It might be intriguing to note that a 2 + 1 dimensional Kondo breakdown quantum critical point (a 
non Fermi liquid) has C ~ T 2 / 3 behavior [30) . Also the same behavior also appear in a class of non- Fermi 
liquids in the large N limit [3Tj . 
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We start with the effective action of Einstein-Maxwell-scalar theory: 
1 



S 



2k 



.2 



-g 



(R-2A)-W((t))F^F f 



fJbV 



l -d^<p-v{4>) 



(5.64) 



-p- is the negative cosmological 



where W(<fi) and V(4>) can be any function. A — ^ 
constant. We would like to mention that this system has been extensively studied in 
[33H35] and some of our solutions, especially our scaling solutions discussed in section 5.1, 
correspond to their special examples. 



5.1 Equation of Motions and General Solutions 

The equations of motion are summarized as follows. The Einstein equation reads 
R^v — -j9\ivR + h-9fiv = 

_EMl F ^F p(T g^ + 2W{(f>)F^F vp - \d p <pd p <Pg, u + \d,<\>d v <\> - ^(0)^.(5.65) 
The Maxwell equation is given by 

d„ {s/^gW{(j))F^) = 0. (5.66) 
Finally, the scalar equation of motion reads 

^(^)_^_»^0. (5.67) 



a/— ~9 d(j) d(j) 

Our ansatz for solutions is given as follows: The metric is again given by f)3.27p . The 
dilaton and the gauge potential are functions of z 

= <f)(z), A t = a{z). (5.68) 

The Maxwell equation ( I5.66P leads to 

A 



a'(z) = W —Jf{z)g{z), (5.69) 

where A is an integration constant. Notice that A is interpreted as the charge density in 
the dual quantum field theory. 

The Einstein equations (15.65!) can be written as 

Gtt = * /w f A+ « + + m^f, 



z 2 J w V 2 7 R 2 W((f)) J v ' Ag(z) 

R 2 . . / V{4>)\ A 2 z 2 . . 1 ,. , 2 



z 2 ^ ' V 2 J R 2 W(<PY y ' 4 
R 2 /. V(0)\ A 2 z 2 1 



G _ = __^ + ^J + ____^. (5.70) 
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The scalar field equation of motion is given by 



■ a , #m aA{z) \ _ ?m + <m+) = p. (5 . 71) 



Now, for our background are found to be 

f(z)(3g{z) + zg'(z)) 



Gu 

G x: 



z 2 g{z) 2 
3f(z)-zf'(z) 
z 2 f(z) ' 

z 2 g{z)f{z) 2 - 4f(zn3g(z) + zg'(z)) + zf{z){zf{z)g'{z) + g(z)(4/'(s) - 2*/" (*))) 

Az 2 f{z) 2 g{z) 2 

(5.72) 



From (I5.72p and (15.701) . the three Einstein equations are equivalently expressed as 
follows 

1/(0) 

*W/'(*) 2 +4/(zfl-6g(*) + 6g(^ 

AR 2 f(z) 2 g{z) 2 

1 R 2 zg{z)f{z) 2 + f(z)(zf'(z)g'(z) + g(z)(Af'(z) - 2zf"(z)) 



W(</>) 8A 2 z 3 f(z) 2 g(z) 2 



kl/ , 2(g(z)f'(z) + f{z)g'{z)) 



<P\zY = - ^ " \ " . (5.73) 

zf{z)g{z) 

It is immediate to find the physical conditions by requiring W(<f>) > and (f)'{z) 2 > 0, 
and we find that they are exactly equivalent to the null energy conditions (I4.53j) 

We can show that (15.711) is automatically satisfied if we substitute (15.731) into it. Thus, 
in principle, we can always have a full solution to all equations of motion by choosing V{4>) 
and W(4>) appropriately. Only non-trivial constraint comes from the no-ghost (and null 
energy) condition (14. 53[) . 

5.2 IR Solution with Logarithmic Entanglement Entropy 

To reproduce (12. 3 j) we need to have 

g(z) = 4, (5-74) 
z F 

in the IR region z » Zp. We also assume the following form of f(z): 

f(z) = kz~ p , (5.75) 
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where k is a positive constant. In this case, the solution exists for p > 2 and it is given by 



00) = y/2(j>-2)]ogz, (5.76) 

T/^ 6 (32 + 12p + p 2 )4 

= ^ ^ , (5-77) 

8A 2 z 6 , . 

w ® = 4p(s + p)R- (5 ' T8) 



Therefore we find 

6 (32 + 12p + p 2 )4 _ 
^ = JP AR* 6 

z 2 F p(p + 8)R 2 



= ..2 7Z , oNm ev/ ^ T ( 5 - 79 ) 



It is clear that this solution gets singular in the IR limit z = oo. In order to show that 
the singularity at zero temperature is not problematic, we need to find regular black hole 
solutions at finite temperature. We assume the potentials V(<p) and W(4>) are given by 
(I5.77P and (I5.78p . Then we can find the following simple black hole solutions 

9(z) 



z 2 F h(z) ' 
f(z) = k-z- p - h(z) 



(z) = v /2(p-2)\ogz, (5.80) 



where h(z) is given by 



p+n 



h(z) = l-Mz—. (5.81) 

The positive parameter M represents the mass of the black hole and the horizon is situated 
at h(r) = 0. In this example, we can show that the specific heat C or equally the thermal 
entropy S behaves as follows 

/ T \ ?T4 

C oc S oc = . (5.82) 



\zpyk J 

This temperature dependence is consistent with the general argument in (I4.60p . 

5.3 Embedding into Asymptotically AdS Solutions 

Finally, we would like to embed the above IR solution into an asymptotically AdS back- 
ground. To make the presentation simple, we concentrate on the p = 3 case below. Then 
we can choose 

k ■ z~ 3 
1 + k- z~ 




g(z) = \\l + — J f(z)= - ; - — . (5.83) 
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Figure 3: The profile of scalar potential V(<p) + 2A for zp = k = 1 and p = 3. The unit 
of the vertical axis is 1/R 2 . The initial condition of the scalar field is set to be 0(0) = 0. 



The scalar field is found to be 



M^l^ A-Z—-! . (5.84) 

In order to keep the inside of the square root positive for all z, we need to require 

k < 2z\. (5.85) 



The scalar field behaves like 



/3k 

4>{z) ~ v^log^ + ^o, (5.86) 

where 0o is a certain constant defined by the above asymptotic behavior. 
The resulting scalar potential V(<f>) is plotted in FigJSJ It behaves like 

^w--^ 2 . (*-><>) < 5 - 87 ) 

w) .|__(^±i|±i!M e -^-« ( ^ oo) . (5 . 88) 

Notice that the tachyonic mass (15.871) marginally satisfies the BF bound. 
The behavior of W is plotted in FigJH It looks like 

W{<fy ~ 2 5 3- n / 3 ^^0- 4 / 3 , (0 -> 0) 
18 
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Figure 4: The profile of W(4>) for zf = k = 1 and p = 3. The unit of the vertical axis is 
A 2 /R 2 . The initial condition of the scalar field is set to be 0(0) = 0. 



As a function z, W is given by 



W{Z) ~ 3^(94 + 2^ + 11^) • ^ 90 > 

It will be an very intriguing future problem to find a string theory embedding which has 
these properties of V and W. 

In summary, in order to have expected solutions, the potentials should behave as 
follows in the limit <p — > oo 

^'-R 2 R 2 
W{<j>) ~ c 2j R 2 e 3 ^, (5.91) 

where the dimensionless constants C\ and c 2 are fixed once we decide a microscopic theory. 
If we set k ~ 0(z F ) motivated by (15.851) . then we find 0o ~ — \/21og zp. Thus we find the 
relation (up to a numerical constant) 

(5.92) 



The gauge potential at the AdS boundary z = is found to be (again up to numerical 
constants) 

poo A p2 i 

- 4,(0) = i dz wir)^^~A4~-, (5 ' 93) 



This relates the Fermi energy \xp = A t (0) to the charge density A by fip oc y A. In this 
case, the specific heat (I5.82|) behaves like 



10 



C(xTr(n F )-. (5.94) 
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